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Abstract. 

We analyze the quality of the quantum information transmission along a correlated 
quantum channel by studying the average fidelity between input and output states 
and the average output purity, giving bounds for the entropy of the channel. Noise 
correlations in the channel are modeled by the coupling of each channel use with 
an element of a one dimensional interacting quantum spin chain. Criticality of the 
environment chain is seen to emerge in the changes of the fidelity and of the purity. 



1. Introduction 

The common scenario in quantum communication protocols is constituted by two 
distant parties, Alice and Bob, who want to exchange information through a quantum 
communication link. Due to unavoidable noise in the channel, this cannot be perfectly 
accomplished, and some strategies aimed to reduce communication errors have to be 
employed. These are based on complex encoding/decoding operations and on suitably 
tailoring the physical system that acts as a channel. In this context the effect of noise 
on the quantum communication is typically quantified by the so called capacities of 
the channel, that is the optimal rates at which (quantum or classical) information can 
be reliably transmitted in the limit of infinite channel uses [1]. The vast majority 
of the results obtained so far focused on the case of memoryless quantum channels, 
where the noise acts independently for each channel use. However, in real physical 
situations, correlations in the noise acting between successive uses can be established. 
When this happens the communication line is said to be a memory channel, or more 
precisely, a correlated channel. The analysis of these setups is much more demanding 
than the memoryless case, and, at present, only a restricted class of them has been 

solved [aSlllElElElElII!. 

Recently, a physical model for representing correlated channels has been proposed 
in Refs. [HllO], which, in the context of Bosonic channels and qubit channels respectively, 
has established a direct connection between these systems and many-body physics. The 
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Figure 1. Model for memory channels: The red dots represent the channel uses 
(ordered, for mstance, starting from the left to the right). The channel noise is 
modeled as a collection of local interactions between the channel uses and the many- 
body environment E (black dots). 



setup discussed in these proposals is depicted in Fig. [T] Here Alice sends her messages 
to Bob by encoding them into a n-long sequence of information carriers S (the red 
dots of the figure), which model subsequent channel uses associated with n independent 
Bosonic modes |1] or n independent spins . The correlated noise of the channel is then 
described by assuming that each carrier interacts independently with a corresponding 
element of a n-party environment E (sketched with the connected black dots of the 
figure), which, in Refs. [1] and [10], represents a multi-mode Gaussian state and a 
many-body spin state, respectively. Thus given an input state ps of the n carriers, the 
corresponding output state associated with the channel is 

Snips) =TTE[Uips0aE)U^], (1) 

where aE is the joint input state of E and the partial trace is performed over the 
environment. In this equation U represents the unitary coupling between S and E, 
which is expressed as 

n 

W = (g) t/W , (2) 

with U^^^ being the interaction between the i-th. carrier and its environmental 
counterpart (in Ref. [Ij these were beam-splitter couplings, while in Ref. [lOj they 
were phase-gate couplings). Within this framework, memoryless channels Sn = 
are obtained for factorizable environmental input states, while correlated noise models 
correspond to correlated environmental states ce- Interestingly enough, in Ref. [10] 
it was shown that it is possible to relate the quantum capacity [TT] of some specific 
channels ([1]) to the properties of the many-body environment E. 

In this paper we discuss a variation of the model ([1]), which allows us to adapt some 
of the techniques used in Ref. pl2j for characterizing the decoherence effects induced by 
spin quantum baths, in order to analyze the efficiency of a class of correlated qubits 
channels. To do so we consider a unitary coupling U that does not factorize as in 
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Figure 2. Generalized model of spin chain memory channels. As an example, in this 
figure we set m — 2. 



Eq. ([2]). Instead we assume to be a spin chain characterized by a free Hamiltonian 
T-Ce, whose elements interact with the carriers S through the local Hamiltonian Timt- 
With this choice we write 

U = exp[-i{nE + nint)t] , (3) 

with the interaction time t being a free parameter of the model. In particular, as a chain 
Hamiltonian, we consider a spin-1/2 XY model in a transverse field, which can exhibit, 
in some parameters region, ground state critical properties that greatly enhance spin 
correlations [13]. Therefore the distance of the chain from criticality is non trivially 
related to memory effects in the channel. In the second part of the paper, we generalize 
the previous scheme by introducing a given number of m extra spins between any two 
consecutive qubits, as shown in Fig. [2l In this case, we can use the number m to 
modulate the memory effects. 

2. The Model 

As the environment E of the system in Fig. [T]we consider an interacting one-dimensional 
quantum spin-1/2 chain described by an XY exchange Hamiltonian in a transverse 
magnetic field: 

n 

= - 2 5Z + ^^^i ^^+1 + - + 2Aa;] , (4) 

where a" (with a = x,y, z) are the Pauli matrices of the j-th spin, J is the coupling 
strength between neighboring spins, and A is the external field strength 0. The model in 
Eq. (jlj) for < 7 < 1 belongs to the Ising universahty class, and has a critical point at 
Ac = 1; for 7 = it reduces to the XX universality class, that is critical for |A| < 1 [13j . 

Following Ref. [12], we then assume that each carrier qubit is coupled to one 
environmental spin element through the coupling Hamiltonian 

^mt(j) = -^^|e)j(e| ®(t| , (5) 
§ Hereafter we always use open boundary conditions, therefore we assume cr"cr" = 0. 
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where \g)j and \e)j respectively represent the ground and the excited state of the j-th 
qubit. Hence the total Hamiltonian H = He + 'Hmt is given by 

J n n 

i=i i=i 



Finally, as in Refs. [TOl [12], we suppose that at time t = the environment chain is 
prepared in the ground state l*/?)^ of T-Le- We then consider a generic input state \ip) ^ 
of the n qubit carriers of the system (i.e., the input state of the red dots in Fig. [T]), and 
write it in the computational basis: 

X 

where are complex probability amplitudes and the sum runs over N = 2^ possible 
choices of x, each of them being a binary string of n elements in which the j-th element 
is represented as g or e, according to the state (ground or excited, respectively) of the 
corresponding j-th qubit. 

For each vector |a;)^ we define as the set of the corresponding excited qubits (for 
instance, given n = 5 and \x) g = \egeeg)g, then Sx contains the 1^*, 3'''^ and 4*^^ qubits). 
After a time t, the global state of the qubits and the chain will then evolve into 

\lp)s^\<f)E ^ '^(^x\x)s®l^x\f)E , (8) 

X 

where U is the global evolution operator of Eq. ([3]), while Ux = exp[— i H% t] is associated 
to the following chain Hamiltonian: 

ni = nE-eY,^h (9) 

According to Eq. (JT]), the channel output state is then described by the density matrix 
^n(l^)5(V^|) = Ps = ^Lxy axa*y \x)s{y\ , (10) 

x,y 

where 

Lxy = e{^\UIUx \'^)e , (11) 

can be seen as a generalized Loschmidt echo, denoting the scalar product of the input 
environment state Iv?)^ evolved with and Uy, respectively [H]. These quantities 
can be evaluated by first mapping the Hamiltonian (Q into a free-fermion model via a 
Jordan Wigner transformation [15] 

Ck = exp f in ^ a^aj j a' , (12) 
V i=i ^ 

where = [a^ ±ia'^)/2, and then by diagonalizing it with a Bogoliubov rotation of the 
Jordan Wigner fermions {c\, Ck}k=i,...,n- This allows one to find an explicit expression 
of the Loschmidt echo in terms of the determinant of a 2n x 2n matrix (see Ref. [12] for 
details): 

Lxy = i^(^|e^^-*e-^^-V)i? = det(I - po + Poe^^'^^e-^^^*) , (13) 
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where Hk = Eij = (4 • • • cjv Ci • • • cn^ , and [po]ii = E{'^\'^\'^j\<p)E are 

the two-point correlation functions of the chain. 



3. The channel 



The echoes ( ITTj) provide a complete characterization of the correlated channel £n- In 
particular, since L^^ = 1 for all x, Eq. ffTOj) shows that the channel £^^1 is unital, i.e. 
it maps the completely mixed state j^^^^ \^)s{x\ into itself. Furthermore, the N x N 
matrix of elements L^y/N coincides with the Choi-Jamiolkowski state [16] of the map. 
The latter is defined as the output density matrix obtained when sending through the 
channel £n half of the canonical maximally entangled state = J2x \^)s ® \^)a 

of the A^-level system S, i.e. 

x,y 

with A being a A^-dimensional ancillary system and X4 being the identity map. Similarly 
to the case analyzed in Ref. [10], this is a maximally correlated state [17] whose 1-way 
distillable entanglement is known to coincides with the "hashing bound" [TTJ [181 IE] ■ 

D^{J{Sn)) = H{Js{£n)) - H{J{£^)) = log2 N - H{J{S^)) , (15) 

where Js{£n) = Tr^[J(£^„)] is the reduced density matrix of J{Sn) associated with 
the system S, while H{-) = Tr[(-) log2(-)] is the von Neumann entropy. At least for 
the subclass of forgetful channels |[8j, the regularized version of Eq. (fT5|) can then be 
used [iHl ED] to bound the quantum capacity [H [TT] of Sn- This is\^ 

n— »oo fl n— +00 n 

The quantity H{J{£n)) corresponds to the entropy of the channel Sn of Ref. [21j, which 
can be used as an estimator of the channel noise. In our case it has also a simple 
interpretation in terms of the properties of the many-body system E: it measures the 
entropy of the ground state! after it has evolved through a random application of 
the perturbed unitaries 151, i.e. 

H{J{Sn)) = H{a'E) , with ct'^^^Y.^, • (1'^) 

X 

II The inequality (|16p is a consequence of the fact that the quantum capacity Q of a channel does not 
increase if we provide the communicating parties with a 1-way (from the sender to the receiver) classical 
side communication line |18[ I20j . It is derived by constructing an explicit quantum communication 
protocol in which i) Alice sends through the channel half of the maximally entangled state |+)5a to 
Bob, ii) the resulting state J{£n) is then 1-way distilled obtaining Z?i(J(£'„)) Bell pairs which, finally, 
in) are employed to teleport Alice messages to Bob. It is worth noticing that for the channel analyzed 
in Ref. [10] the right hand side of Eq. (|15p was also an upper bound for Q. 

^ This is a trivial consequence of the fact that J{£n) is the reduced density matrix of the pure state 
Tlx \^^)sA ®Ux\^p) E tracing out the environment, and of the fact that the von Neumann entropies 
of the reduced density matrices of a pure bipartite system coincide. 
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Figure 3. Averaged channel fidelity as a function of the interaction time, for different 
values of the transverse field A: from right to left A = 0.25, 0.5, 0.75, 0.9, 1, 1.1, 1.5, 2. 
Here we simulated a channel of n — 50 qubits coupled to an Ising chain, and set 
an interaction strength e = 0.05; the fidelity has been evaluated by sampling over 
Nav = 5 X 10* randomly chosen initial conditions. In the inset we plot the averaged 
fidelity at a fixed interaction time Jt* = 5, as a function of A. 



Unfortunately, for large n the computation of the von Neumann entropy of the state 
J{Sn) is impractical both analytically and numerically, since it requires to evaluate an 
exponential number of L^y elements. Interestingly enough, however, we can simplify our 
analysis by considering the fidelity between J{Sn) and its input counterpart |+)5a (see 
Eq. ( fT4l) ). As discussed in the following section, this is a relevant information theoretical 
quantity, since it is directly related to the average fidelity between input and output 
state of the channel £n and provides us an upper bound for H{J{Sn))- Similarly we can 
compute the purity of J{Sn) which, on one hand, gives a lower bound for H{J{£n)), 
while, on the other hand, it is directly related to the average channel output purity of 
the map £n- 

4. Average transmission fidelity 

According to Eq. f|T^ . the fidelity between the Choi-Jamiolkowski state J{Sn) and its 
input counterpart coincides with the average value of the Loschmidt echoes L^-y, i.e. 
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Even without computing all the L^y, this quantity can be numerically evaluated by 
performing a sampling over A^^^ randomly chosen couples [x, y) of initial conditions, 
and averaging over them 

J^^J^^^= J2 MLxy], (19) 

{x,y)=l 

(where we used the fact that Ly^ = L*^y)- The quantity JF provides an upper bound for 
H{J{Sn)) through the quantum Fano inequahty [22], i.e. 

H{J{S„)) ^ H^iJ") + (1 - -F) log2(4" - 1) ^ H^iJ') + 2n{l - ^) , (20) 

where i^2(') = ~(")log2(') ~ ~ (■)]log2[l ~ (')] is the binary entropy functior0. 
Furthermore is directly related to the average transmission fidelity {F) of the map 
En- For a given pure input state ([7]), the transmission fidelity is 

F(v^) ^ smn^sm \i^)s = 5^i:x,i«xiv,r . (21) 

x,y 

Taking the average with respect to all possible inputs, we get 

{F)=J2LxyPxy, (22) 

x,y 

where pxy = with (...) being the average with respect to the uniform Haar 

measure. The probability distribution p^y can be computed by using simple geometrical 

1+2 5x,y 

N{N+2) 



arguments [23]. As shown in Appendix A this yields pxy = ]^('^l2) and hence: 



where we used the fact that L^x = 1 and L^y = L*^. Therefore, from Eq. (ITSll we get 

The fidelities JF and (F) are not directly related to the channel quantum capacity, 
nonetheless, as in Eq. fl20|) . they can be used to derive bounds for Q % More 
generally, values near to unity of the fidelity between the output channel states and 

We numerically checked the convergence of J- with Nav We first considered a situation with a 
few number of qubits {n < 10), such to compare sampled averages, J-av, with exact averages over all 
possible events, Tex — J' ■ We found that, already at N^v = 10*, absolute differences \Tav — Tex\ are 
always less than 2 x 10~^, while at Nav = 5 x 10'' the error is less than 5 x 10"'^, independently of 
the values of the interaction time t, the transverse field A and the system size n. In a second time, we 
simulated systems with definitely larger sizes (n « 50) and simply check the convergence of Tav with 
Nav Differences between fidelities with Nav = 10* and Nav = 5 x 10* are of the same order as the 
deviation of the curve with Nav — 10* from the exact one for small sizes. Therefore we can reliably 
affirm that fidelity results with Nav = 5 x 10* are exact, up to an absolute error of order 5 x 10~^. 
* Equation ([20|l can be easily derived by noticing that H{J{£n)) and T coincide, respectively, with 
the exchange entropy and entanglement fidelity of the channel £n associated with the maximally mixed 
state Is/2" of S". 

jj In particular from Eq. (|20p and Eq. (fT6|l one gets Q ^ 1 — 2 lim„^oo J^. 




Figure 4. Fidelity for a channel coupled to an Ising chain with A = 1 and different 
qubit numbers n: from right to left n = 4, 6, 8, 10, 16, 30, 50; the interaction strength 
is kept fixed at e = 0.05; data are averaged over Nav = 5 x 10"* configurations. Inset: 
Tav at a fixed time t*, as a function of n. 



their corresponding input states, are indicative of a fairly noiseless communication line. 
On the contrary, values of the transmission fidelities close to zero, while indicating 
output states nearly orthogonal to their input counterparts, do not necessarily imply 
null or low capacities, since such huge discrepancies between inputs and outputs could 
still be corrected by a proper encoding and decoding strategy (e.g. consider the case of 
a channel which simply rotate the system states). 

Equation (|T3i) allows us to numerically compute the averaged transmission 
fidelity f|T9l) : numerical results in this and in the next sections are given for the case 
7 = 1, i.e. we study the correlated quantum channel defined by the Ising model. The 
behavior of the averaged channel fidelity JF^^,, defined in Eq. (fT9|) and related to the 
memory channel scheme of Fig. [H with respect to the interaction time t (the free model 
parameter) is shown in Fig. [3j the plots are given for n = 50 qubits, each of them 
coupled to one spin of an Ising chain with a coupling strength e = 0.05. Different 
curves stand for different values of the transverse magnetic field A: as it can be clearly 
seen, the fidelity J^av decays as a Gaussian in time, irrespective of the field strength 
A 0. The signature of criticality in the environment chain can be identified by studying 
the function jFat,(A,t) for fixed interaction time t*: the inset of Fig. [3] displays a non 

fl For times longer than those in the scales of Figs. [3] and[4j revivals of the fidelity are present. See 
[Appendix B| for details. 
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Figure 5. Short-time Gaussian decay rate a as a function of the transverse field A. 
We rescaled a with the system size according to the scahng of Fig. HI In the inset we 
plot the first derivative of the same curves in the main panel, with respect to A. 

analytic behavior for the derivative of J-'av{^,t*) with respect to A at the critical point 
Ac = 1. This non analyticity will be clearer in the following, where the size scaling will 
be considered. In Fig. H] we study the behavior of J-'av with the number n of qubits, 
at a fixed value of transverse magnetic field A = 1. As it is shown in the inset, at a 
given interaction time t*, the fidelity J-'av{^c,t*) decays exponentially with n (the same 
behavior is found when A 7^ 1). This dependence of the decay rate implies that the 
average fidelity J-'av can be fitted by: 

J^av ~ e~°*^ with a oc n ; (25) 

in other words, the Gaussian decay rate is extensive. Indeed, since the fidelity is a global 
quantity that describes the evolution of the state of the whole n-body system, it should 
start decaying as a Gaussian (at least at small times) [21], with a reasonably extensive 
decay ratio. This prediction is confirmed by the results of Fig. [5l where we report the 
decay rate as a function of the transverse magnetic field for different system sizes n. In 
proximity of the critical point, the decay rate undergoes a sudden change, which becomes 
more evident when increasing the system size. The signature of criticality at Ac = 1 and 
the finite size effects can be better analyzed by looking at the derivative of the decay rate 
with respect to the transverse field: The inset of Fig. [5] clearly show that dxa exhibits 
a non analytic behavior at the critical point Ac, at the thermodynamical limit. Notice 
also that, due to finite size effects, the maximum of dxa does not coincide exactly with 
the critical point, that can be rigorously defined only at the thermodynamical limit, but 
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occurs at a slightly smaller value of A. However, we checked that a finite size scaling 
gives the right prediction of the critical point located at A = 1. 

5. Average output purity 

Another quantity that can be evaluated with relatively little numerical effort is the 
purity of the Choi-Jamiolkowski state J(£^„), i.e. 

P,^Tr[J(^„)2] = ^5^|L.,r (26) 

xy 

As in the case of JF, this can be computed by approximating the summation with a 
random sampling, i.e. 

V2^Va.^— \^-y\'- (27) 

(x,y)=l 

The quantity (l26l) provides us two important pieces of information. First of all, it yields 
a useful bound on the channel entropy H{J{Sn))- This follows from the inequality [25] 

H{J{Sn)) ^ H2{J{Sn)) = -\og^V2 , (28) 

with H2{-) = — log2 Tr[(-)^] being the Renyi entropy of order 2. Furthermore V2 is 
directly related to the average output purity (P2) of the channel £n- This is obtained by 
averaging over all possible inputs \4')s the purity of the output state SnU'ip) sii'l), i-e. 

(A) = (Tr [{SMsm?] ) = E \^^y\" (l"^l>.l') = E \^^y\"P^y ' (29) 

x,y x,y 

where we used Eq. ( JTOl) and where pxy are the probabilities defined in Eq. ( ]A.2l) . 
According to Eq. (12^ this yields, 

<^^> = wrf^- + nt-2 ■ <^°' 

The average purity is a rather fair indicator of the noise induced by the coupling to the 
environment: if the carrier qubits get strongly entangled with the environment, P2 is 
greatly reduced from the unit value (for large n it will tend to zero); on the other hand, 
a channel which simply unitarily rotates the carrier states has a unit purity. However, 
we should stress that also the purity may intrinsically fail as a transmission quality 
quantifier: there are strongly noisy channels with very high output purity (consider, for 
example, the channel which maps each input state into the same pure output state). 

Then, we study the average channel purity Vav as a function of the model free 
parameter, the interaction time t. The results are reported in Fig. [6]for a chain of n = 30 
qubits and different values of the transverse field A. We notice qualitatively different 
behaviors depending on the values of the transverse field A: If A < Ac the averaged 
purity oscillates in time and asymptotically tends to an average constant value; as far 
as the critical point is approached, Vav drops to smaller values (revivals are again due to 
finite system size effects), reflecting the fact that at criticality correlations between the 



Spin chain model for correlated quantum channels 11 




Figure 6. Averaged purity of tlie channel output state as a function of the interaction 
time, for various transverse field strengths: from left to right A = 0.25 (black circles), 
0.5 (red squares), 0.75 (green diamonds), 0.9 (brown triangles up), 1 (blue triangles 
down), 1.1 (magenta triangles left), 1.5 (orange triangles right), 2 (violet crosses). Here 
we simulated a channel of n = 30 qubits coupled to an Ising chain, set an interaction 
strength e = 0.05, and averaged over Nav = 5 x 10'* random initial conditions. In the 
inset we plot the short-time Gaussian decay rate /? as a function of A. 

qubits and the environment are stronger. Crossing the critical point, in the ^-ordered 
phase (A > Ac) the purity is generally higher and asymptotically takes values very close 
to the unit value. This can be easily understood in the limit A — > cxd: in this case the 
spins in the chain are "freezed" along the field direction and they cannot couple with 
anything else, resulting in a watchdog-like effect [2^. Independently of the transverse 
field value, the average purity Vav decays as a Gaussian Vav ~ e~^*^ in the short time 
limit. As for the averaged fidelity J-'av, we have then analyzed the decay rate /5 as a 
function of A: as before, f3 exhibits a signature of criticality via a divergence, at the 
thermodynamical limit, in its first derivative with respect to A (see the inset of Fig. [6]). 

6. Generalized model 

We finally concentrate on the generalized model depicted in Fig. [21 where a certain 
number m of environment spins are present between two consecutive spins coupled to 
the qubits. The richness of the model, that is characterized by a large number of 
parameters, and by a global size which grows both with n and m, requires a huge 
numerical effort in order to simulate it, therefore we decided to analyze only the average 
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Figure 7. Upper panel: average channel fidelity for the generalized model in Fig. [51 
with n = 12 qubits, e = 0.05 (averages have been performed over Nav = 
initial states). The various curves are for different numbers m of spins between two 
consecutive qubits, and different values of transverse magnetic fields A = 0.25, 1, 2. 
Lower panel: absolute differences in the fidelities between configurations at various m, 
as a function of A and at a fixed interaction time Jt* = 10. 



channel fidelity J-'av In the upper panel of Fig. [7] we show J^av as a function of the 
interaction time t for different values of m and three values of the transverse field 
A = 0.25, 1, 2; we fix a number of qubits n = 12 and an interaction strength e = 0.05. 
Hereafter we will concentrate on this typical result, as we performed some 

checks with larger numbers of channel uses (n = 30, 50 and m = 0, 1, 2), and found 
qualitatively analogous results. We immediately observe that differences for various m 
are tiny, even if the fidelity generally tends to increase when increasing m; the sensibility 
with m suddenly enhances at criticality (Ac = 1), where correlations in the environment 
decay much slower than in the other cases. On the other hand, when A is far from 
Ac, differences between fidelities upon a variation of m are greatly suppressed, and the 
generalized model mostly behaves as the model in Fig. [H Again, this reflects the fact 
that, out of criticality, each qubit is mostly influenced only by the spin that is coupled 
to, as the spin does not exchange correlations with the other environmental spins. The 
resulting channel properties are then defined only by the local properties of the chain. 
On the contrary, at criticality, the spins are correlated and then the resulting channel 
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Figure 8. Gaussian decay rate of tlie averaged fidelity for model in Fig. [2l as a 
function of the transverse field A. The various curves stand for different numbers of 
spins m between two consecutive qubits (here n = 12, e = 0.05, Nav = 10^). In the 
inset we show the first derivative with respect to A of the curves in the main panel. 



properties are influenced by the distance of the spins coupled with the qubits. In the 
lower panel of Fig. [7] we explicitly plot the differences in the fidelities for various m as a 
function of A, and for a fixed interaction time; a peak in proximity of Ac is clearly visible. 
We point out that, as already noted at the end of Sec. H] concerning the size scaling of 
the fidelity, the maximum in the differences does not occur exactly at the critical point. 

The sensitivity to criticality is again demonstrated by the averaged fidelity Gaussian 
decay rate a as a function of A, as shown in Fig. [8] for different values of m: the first 
derivative in the inset has a maximum in correspondence of a value that approaches 
the critical point Ac at the thermodynamical limit. Indeed, increasing m is equivalent 
to approaching the thermodynamical limit of the chain, thus resulting in an increase 
of the quantum phase transition effects. A double check of this comes from the cyan 
triangles- down curve of Fig. (8) in this case we take m = 4, but we break one of the 
links between two intermediate spins. The environment is then formed by disconnected 
chains, each of them made up by 5 spins, therefore the system cannot undergo a phase 
transition in the limit n oo: the signature of criticality has completely disappeared. 

7. Conclusions 

In conclusion we have introduced and characterized a class of correlated quantum 
channels, and we have given bounds for its entropy by means of the averaged channel 
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fidelity J^av and purity Vav Even though in general these bounds might not be strict, we 
give a characterization of the channel in terms of quantities that have a clear meaning 
from the point of view of the many body model we have introduced. 

In the case of an environment defined by a quantum Ising chain, we have shown that 
the averaged channel purity and the fidelity depend on the environment parameters and 
are strongly influenced by spin correlations inside it, in particular by the fact whether 
the environment is critical or not. We expect that some different environment models, 
such as, for example, the XY spin chain, will behave qualitatively similarly of what 
found in this work, as it belongs to the same universality class. This might not be the 
case for other models, like the Heisenberg chain, which will be object of further study 
in the near future. 
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Appendix A. 

The probability p^y of Eq. (122|) can be computed as follows: we first define = 
and convert the string x into a decimal number from 1 to = 2" by trivially identifying 
g = and e = 1. The average over a uniform distribution of all pure input states on 
the Bloch hypersphere for n qubits is 

P.y = CN [ dn---/ drNrlrl6{l-r^), (A.l) 
Jo Jo 

where = rf + ■ ■ ■ + r^, and C^^ = dri ■ ■ ■ Jq drN 6{1 — r^) is a normalization 
constant. Changing the limits of integration due to the delta function and using 
the property of the Gamma function T{z) = f^"^ y^^^e^'^dy, it is easy to show that 
Cn = 2^7r-^/2r(Ar/2) and 

p = ^ + '^^-^y . (A.2) 
^"^y N{N + 2) ^ ^ 



Appendix B. 

For times longer than those in the scales of Figs. [3] and HI time revivals of the fidelity 
are present, i.e. the fidelity increases back towards the unit value periodically, due to 
the finite system size. For finite values of the transverse field revivals are not perfect, 
that is J-'av{t) 7^ 1 for t > 0. Anyway, as far as A increases, the revivals are stronger 
and happen with period t^ which does not depends on the system size n. This can be 
understood in the limit A +oo, where the ground state of the environment If)^ is a 
fully ^-polarized state, thus being an eigenstate of the chain Hamiltonian T-C% in Eq. ([9]): 
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the generalized Loschmidt echo of Eq. (ITT]) is then given by L^y = e"''^-'^^-^'*, where N'x-y 
is the number of excited qubits in the sequence x minus the one in the sequence y. It 
is easy to see that there are 2" x (2) different possibihties to choose two sequences x, y 
such that the corresponding states differ in the state of k qubits, then Afx-y = ±i with 
j = 0, 2, . . . , A; (if is even) or j = 1,3, ... ,k (if is odd). Therefore, when averaging 
over input states, each term contributes with Pxy^^^^^^- Noting that (U^UxY = UlUy, 
we have 

(F) = 2'^pxx + Pxy (co + Ci cos(£:t) + C2 cos(2£:t) + . . . + c„ cos{net)) , (B.l) 
It follows a perfect revival for the fidelity at times t^ such that et^ = 2tt [mod 2tt]. 
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